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Bacillus  anthracis  is  the  etiologic  agent  of  anthrax,  an  acute  fatal 
disease  among  mammals.  It  was  thought  to  differ  from  Bacillus 
cereus,  an  opportunistic  pathogen  and  cause  of  food  poisoning,  by 
the  presence  of  plasmids  pXOI  and  pX02,  which  encode  the  lethal 
toxin  complex  and  the  poly-y-D-glutamic  acid  capsule,  respectively. 
This  work  describes  a  non-B.  anthracis  isolate  that  possesses  the 
anthrax  toxin  genes  and  is  capable  of  causing  a  severe  inhalation 
anthrax-like  illness.  Although  initial  phenotypic  and  16S  rRNA 
analysis  identified  this  isolate  as  B.  cereus,  the  rapid  generation 
and  analysis  of  a  high-coverage  draft  genome  sequence  revealed 
the  presence  of  a  circular  plasmid,  named  pBCXOI,  with  99.6% 
similarity  with  the  B.  anthracis  toxin-encoding  plasmid,  pXOI. 
Although  homologues  of  the  pX02  encoded  capsule  genes  were 
not  found,  a  polysaccharide  capsule  cluster  is  encoded  on  a  second, 
previously  unidentified  plasmid,  pBC218.  A/J  mice  challenged  with 
B.  cereus  G9241  confirmed  the  virulence  of  this  strain.  These 
findings  represent  an  example  of  how  genomics  could  rapidly 
assist  public  health  experts  responding  not  only  to  clearly  identi¬ 
fied  select  agents  but  also  to  novel  agents  with  similar  pathogenic 
potentials.  In  this  study,  we  combined  a  public  health  approach 
with  genome  analysis  to  provide  insight  into  the  correlation  of 
phenotypic  characteristics  and  their  genetic  basis. 

Bacillus  cereus  and  Bacillus  anthracis  are  members  of  a  closely 
related  phylogenetic  cluster  referred  to  as  the  B.  cereus 
group.  Although  classification  of  B.  anthracis  and  B.  cereus  as 
separate  species  based  on  molecular  analysis  has  been  ques¬ 
tioned,  they  differ  in  phenotype  and  in  the  diseases  they  cause 
(1-4).  B.  cereus  is  ubiquitous  in  nature  and  is  an  opportunistic 
pathogen.  In  individuals  with  underlying  conditions,  immuno¬ 
compromised  individuals,  or  patients  recovering  from  surgery, 
B.  cereus  has  been  known  to  cause  a  variety  of  infections, 
including:  endophthalmitis,  bacteremia,  septicemia,  endocardi¬ 
tis,  salpingitis,  cutaneous  infections,  pneumonia,  and  meningitis 
(5).  In  contrast,  B.  anthracis  is  the  etiologic  agent  of  anthrax,  an 
acute  fatal  disease  among  mammals,  which  in  recent  years  has 
become  known  for  its  use  as  a  biological  weapon  (6).  The  genes 
encoding  the  major  anthrax  toxins  and  the  poly-y-D-glutamic 
acid  capsule  are  located  on  two  virulence  plasmids,  pXOI  (182 
kb)  and  pX02  (96  kb),  respectively,  and  are  required  for  full 
virulence  (7-9). 

Materials  and  Methods 

Genome  Sequencing.  Genomic  DNA  extraction  of  B.  cereus  G9241 
was  performed  at  the  Centers  for  Disease  Control  and  Preven¬ 
tion.  Random  shotgun  sequencing  was  undertaken  by  using  a 
strategy  similar  to  previous  projects  (10).  Random  insert  librar¬ 
ies  of  2. 5-4.5  kb  and  4.5-9. 0  kb  were  constructed,  and  93,562 


high-quality  sequences  were  obtained  (817-bp  average  read 
length).  A  13.5 X  coverage  draft  genome  sequence  was  obtained 
and  assembled  by  using  the  Celera  assembler  (11).  The  assembly 
consisted  of  131  contigs  that  were  ordered  and  concatenated 
with  our  prior  knowledge  of  the  B.  anthracis  and  B.  cereus 
genome  sequences  (10, 12).  The  resulting  5,286,464-bp  sequence 
was  annotated  through  the  Institute  for  Genomic  Research 
Bioinformatics  pipeline  as  reported  (13).  The  genomic  features 
of  this  molecule  are  described  in  Table  2,  which  is  published  as 
supporting  information  on  the  PNAS  web  site. 

Capsule  Visualization.  B.  cereus  was  streaked  onto  trypticase  soy 
agar  containing  5%  sheep  blood  (Becton  Dickinson  Microbiol¬ 
ogy  Systems)  and  incubated  in  ambient  atmosphere  overnight  at 
37°C.  B.  anthracis  was  streaked  onto  trypticase  soy  agar  con¬ 
taining  0.8%  sodium  bicarbonate  and  incubated  in  5%  CO2  at 
37°C.  Cells  from  a  single  colony  were  dispersed  into  5  pi  of  water 
on  a  microscope  slide  and  covered  with  a  coverslip.  A  single  drop 
of  india  ink  (Remel,  Lenexa,  KS)  was  applied  to  the  edge  of  the 
coverslip  and  allowed  to  diffuse  into  the  sample.  The  cells  were 
visualized  with  a  X100  oil-immersion  objective  (14). 

Biochemical  and  Microbiological  Phenotypic  Characterization.  Motil¬ 
ity  was  observed  by  wet  mount  of  cells  grown  in  heart  infusion 
broth  and  by  using  Motility  Test  medium  with  triphenyltetrazo- 
lium  chloride  (Remel).  Susceptibility  to  y-phage  (15)  was  per¬ 
formed  by  adding  5  pi  of  y-phage  to  the  first  and  second 
quadrants  of  a  5%  sheep  blood  plate  streaked  for  isolation  and 
incubated  at  37°C  overnight.  The  B.  anthracis  specific  two- 
component  direct  fluorescent  Ab  (DFA)  assay  was  performed  as 
described  by  De  et  al.  (16).  Sequencing  of  the  16S  rRNA  gene 
was  performed  as  described  by  Sacchi  et  al.  (17). 

Comparative  Genomic  Hybridization  (CGH).  CGH  was  performed  as 
described  for  other  B.  cereus  isolates  by  Read  et  al.  (10). 

Multiple  Locus  Sequence  Typing.  Phylogenetic  tree  of  members  of 
the  B.  cereus  sensu  lato  group  was  generated  based  on  partial 
nucleotide  sequences  of  seven  housekeeping  genes,  totaling 


Abbreviations:  CGH,  comparative  genomic  hybridization;  PA,  protective  antigen;  DFA, 
direct  fluorescent  Ab. 

Data  deposition:  The  Whole  Genome  Shotgun  project  reported  in  this  work  has  been 
deposited  in  the  GenBank  database  (accession  no.  AAEK00000000).  The  version  described 
in  this  paper  is  the  first  version,  AAEK01 000000. 
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2,977  nucleotides  and  was  built  by  using  the  neighbor-joining 
method  (18)  applied  to  a  matrix  of  pairwise  differences  among 
sequences.  Detailed  multiple  locus  sequence  typing  methods  are 
available  online  at  http://pubmlst.org/bcereus. 

blast  Score  Ratio  Analysis.  The  blast  score  ratio  analysis  is  a 
modification  of  the  technique  described  by  Read  et  al.  (19).  For 
each  of  the  predicted  proteins  of  B.  cereus  G9241,  we  obtained 
a  blastp  raw  score  for  the  alignment  against  itself 
(REF_SCORE)  and  the  most  similar  protein  (QUE.SCORE)  in 
each  of  the  genomes  of  B.  cereus  ATCC10987  (12),  B.  cereus 
ATCC14579  (20),  andR.  anthracis  Ames  (10).  These  scores  were 
normalized  by  dividing  the  QUE_SCORE  obtained  for  each 
query  genome  protein  by  the  REF_SCORE.  Proteins  with  a 
normalized  ratio  of  <0.4  were  considered  to  be  nonhomologous. 
The  normalized  blast  score  ratio  of  0.4  is  generally  similar  to 
two  proteins  being  «=30%  identical  over  their  entire  length. 

blast  score  ratios  were  plotted  as  x,y  coordinates  as  shown  in 
Fig.  5,  which  is  published  as  supporting  information  on  the  PNAS 
web  site.  Each  protein  in  the  reference  genome  ( B .  cereus  G9241) 
was  grouped  according  to  its  scores  in  each  of  the  query  genomes 
and  colored  as  follows:  yellow,  unique  to  the  reference;  red, 
common  to  all  three;  cyan,  common  between  query  1  and  the 
reference  but  absent  in  query  2;  blue,  common  between  query  2 
and  the  reference  but  absent  in  query  1. 

Fig.  1C  shows  a  synteny  plot  in  which  each  triangle  indicates 
a  single  protein  plotted  on  the  5'  end  of  the  coding  region  in  the 
reference  genome  and  the  best  match  in  the  query  organism. 

Protective  Antigen  (PA)  ELISA.  Rows  of  an  ELISA  plate  were 
coated  overnight  in  PBS  with  B.  anthracis  PA-specific  mAb 
(Biological  Defense  Research  Directorate)  or  normal  mouse 
IgG  (Sigma)  (negative  capture).  Plates  were  blocked  for  1  h  with 
5%  skim  milk  at  37°C.  After  washing  three  times  with  wash 
buffer  (IX  PBS  with  0.1%  Tween  20),  the  antigen  (i.e.,  super¬ 
natant  of  cells  grown  in  horse  serum)  was  added  to  well  1  and 
titered  across  the  plate  from  wells  1-9  with  1:2  serial  dilution  in 
5%  skim  milk  with  0.1%  Tween  (dilution  buffer).  Plates  were 
incubated  for  1  h  at  37°C  and  then  washed  three  times  with  wash 
buffer.  The  biotinylated  detector  Ab  was  then  applied  in  the 
dilution  buffer.  The  detector  for  this  assay  is  the  same  mAb  as 
the  capture,  only  biotinylated.  The  plate  was  again  incubated  for 
1  h  at  37°C,  then  washed  three  times.  A  streptavidin/horseradish 
peroxidase  conjugate  in  dilution  buffer  was  applied,  and,  after 
incubation  for  1  h,  the  plate  was  washed  three  times.  Substrate 
[2,2'-azinobis(3-ethylbenzthiazoline-6-sulfonic  acid)]  was 
added,  and  the  plate  was  incubated  for  30  min  at  37°C,  at  which 
point  it  was  read  in  an  Ultra  Evolution  ELISA  plate  reader 
(Tecan,  Durham,  NC).  The  positive  cutoff  was  calculated  to  be 
three  times  the  standard  deviation  of  the  mean  of  the  negative 
capture  plus  0.150.  Purified  B.  anthracis  PA  protein  (2-500 
ng/ml;  Biological  Defense  Research  Directorate)  was  used  as  a 
positive  control  and  to  generate  the  standard  curve  for  quanti¬ 
tation. 

Mice  i.p.  Challenge.  Spores  of  B.  anthracis  Sterne  strain,  B.  cereus 
ATCC10987,  and  B.  cereus  G9421  were  prepared  by  the  out¬ 
growth  on  sporulation  medium  (per  liter,  23  g  of  nutrient 
agar/0.025  g  MNSO4-H2O/0.25  g  KH2PO4).  Spore  titer  (and 


Fig.  1.  ( A-B )  Circular  representations  of  the  plasmids  of  B.  cereus  G9241  .The 

plasmids  were  compared  by  using  blastp  to  the  B.  anthracis  pXOI  and  pX02 
plasmids  and  B.  cereus  ATCC10987  pBC10987  and  Bacillus  thuringiensis  israe- 
liensis  pBtoxis.  The  inner  two  rings  represent  the  plasmid  coding  regions  in 
forward  and  reverse  orientation;  color  scheme  is  as  described  (1 3).  The  outer¬ 
ring  color  scheme  is  as  follows:  red,  B.  anthracis  plasmids;  blue,  B.  cereus 
ATCC10987  plasmid  pBC10987;  cyan,  proteins  equally  similar  in  both  species; 
and  green,  proteins  unique  to  that  plasmid.  The  third  ring  represents  the 


remaining  physical  gaps  of  the  plasmids.  (O  Synteny  plots  of  the  plasmids  of 
B.  cereus  G9241  pBCXOI  versus  B.  anthracis  Ames  pXOI.  Each  triangle  rep¬ 
resents  a  single  protein  of  G9241.  They  are  color  coded  as  follows:  red, 
normalized  score  >0.95;  orange,  normalized  scores  >0.80  and  <0.95;  yellow, 
normalized  scores  >0.60  and  <0.80;  green,  normalized  scores  >0.40  and 
<0.60;  cyan,  normalized  scores  >0.20  and  <0.40;  and  blue,  normalized  scores 
>0.0  and  <0.20. 
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percentage  of  spores  to  vegetative  cells)  was  determined  by 
direct  counting  under  a  phase-contrast  microscope.  Spore  sus¬ 
pensions  were  washed  and  diluted  appropriately  in  PBS.  Em¬ 
ploying  the  method  outlined  by  Beedham  et  al.  (21),  six  groups 
of  10,  6-  to  8-week-old,  female,  A/J  mice  weighing  20-24  g 
received  i.p.  injections  either  a  low  (range  1.1  X  104  to  2.7  X  104 
colony-forming  units)  or  high  (range  1.1  X  106  to  5.3  X  106 
colony-forming  units)  dose  of  spores  of  either  the  Sterne  strain 
of  B.  anthracis,  B.  cereus  G9421,  or  B.  cereus  ATCC10987.  The 
A/J  strain  was  selected  because  it  was  found  to  be  sensitive  to 
the  vaccine  strain  B.  anthracis  Sterne  (21).  Spores  for  injection 
were  suspended  in  100  pi  of  PBS.  The  control  group  was  injected 
with  sterile  PBS.  Animals  were  monitored  four  times  daily  for  14 
days  for  signs  of  morbidity/mortality.  Mice  found  unresponsive, 
paralyzed/unresponsive,  or  moribund  were  killed  with  CO2,  and 
blood  and  tissue  (heart,  lung,  liver,  and  spleen)  samples  were 
obtained.  For  mice  found  dead,  only  tissue  samples  were  ob¬ 
tained.  At  day  14,  all  remaining  mice  were  killed  with  CO2  and 
processed  as  described  for  sick  mice  above.  Either  blood  films  or 
touch  preparations  from  heart  tissue  were  performed  for  stain¬ 
ing.  Heart  tissue  and  the  blood  removed  from  the  animals 
postmortem  was  plated  on  5%  (vol/vol)  sheep  blood  agar  and 
incubated  at  37°C  for  24  h  to  confirm  bacterial  infection  as  cause 
of  death.  Student’s  t  test  was  applied  to  determine  the  signifi¬ 
cance  of  the  average  survival  between  two  groups.  This  study  was 
carried  out  in  strict  accordance  with  the  Animals  (Scientific 
Procedures)  Act  of  1986. 

Results 

Isolation  of  B.  cereus  G9241.  During  a  retrospective  analysis  of 
Bacillus  isolates  that  caused  unusually  severe  disease,  we  iden¬ 
tified  an  isolate  (G9241)  that  contained  a  homologue  of  the  B. 
anthracis  pXOl-encoded  PA  gene,  pagA,  but  not  the  pX02- 
encoded  poly-y-D-glutamic  acid  capsule  biosynthetic  genes. 
G9241  was  isolated  from  the  sputum  and  blood  of  a  patient  with 
life-threatening  pneumonia,  whose  history,  clinical  features,  and 
laboratory  findings  resemble  those  reported  for  10  bioterrorism- 
associated  inhalation  anthrax  patients  in  2001  (23).  The  patient, 
like  the  10  inhalation  anthrax  patients  from  2001  (23),  was  not 
immunocompromised,  did  not  take  drugs,  was  a  nonsmoker,  and 
had  no  predominant  underlying  disease.  In  addition  to  fever, 
chills,  and  hemoptysis,  the  patient  sought  care  only  2  days  from 
onset  of  symptoms  (3.5  days  for  2001  patients),  reported  nausea 
and  vomiting  (80%  of  2001  patients),  and  had  an  abnormal  chest 
x-ray  (multiple  abnormalities  in  all  2001  patients).  His  initial 
white  blood  cell  count  (WBC)  was  12,000  per  mm3  (median 
WBC  9,900  per  mm3  for  2001  patients),  and  peak  WBC  was 
22,400  per  mm3  (median  peak  WBC  26,400  per  mm3  for  2001 
patients).  The  patient  received  intensive  care,  mechanical  ven¬ 
tilation  for  44  days,  treatment  with  numerous  antimicrobial 
drugs  (timentin,  erythromycin,  rifampin,  isoniazid,  and  bactrim) 
and  corticosteroids,  and  a  partial  lobectomy  before  he  was 
discharged  and  subsequently  recovered  (Table  3,  which  is  pub¬ 
lished  as  supporting  information  on  the  PNAS  web  site). 

Typing  of  B.  cereus  G9241.  This  anthrax-like  clinical  presentation 
was  attributed  to  an  isolate  (G9241)  identified  as  B.  cereus  by 
using  traditional  phenotypic  characteristics  and  biochemical 
reactions  at  the  Centers  for  Disease  Control  and  Prevention 
Special  Bacteriology  Reference  Laboratory  (Table  1).  By  using 
the  16S  rRNA  typing  scheme  of  Sacchi  et  al.  (17),  it  was 
determined  that  G9241  was  type  7  (GenBank  accession  no. 
AY425946),  one  of  six  16S  types  exclusively  associated  with  B. 
cereus  isolates.  Additionally,  strain  G9241  was  included  in  a 
multilocus  sequence  typing  scheme  that  used  seven  house¬ 
keeping  genes  and  has  been  developed  for  the  B.  cereus  group 
(http://pubmlst.org/bcereus).  This  analysis  showed  that  strain 
G9241  is  closely  related  to  but  distinct  fromB.  anthracis  (Fig.  2). 


Table  1.  Biochemical  and  microbiological  characteristics  of  B. 
anthracis,  B.  cereus  spp.,  and  B.  cereus  G9241 
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*B.  cereus  species  have  a  variable  reaction  to  the  6.  anthracis  cell-wall  DFAtest. 
^Typing  for  16S  rRNA  is  based  on  the  system  of  Sacchi  et  al.  (17).  To  note,  B. 
cereus  spp.  have  a  variety  of  types;  however,  type  6  has  never  been  identified 
in  any  group  other  than  B.  anthracis. 

^B.  anthracis  produces  a  poly-y-D-glutamic  acid  capsule,  and  a  putative  capsule 
operon  has  been  identified  in  B.  cereus  ATCC10987  (12)  and  6.  cereus 
ATCC14579  (20).  However,  under  the  experimental  conditions  employed  in 
this  study,  no  capsule  has  been  observed  in  association  with  either  strain.  6. 
cereus  G9241  expresses  a  capsule  (Fig.  3)  that  is  not  a  poly-y-D-glutamic  acid 
capsule. 

This  conclusion  was  also  supported  by  the  CGH,  which  also 
revealed  a  high  degree  of  hybridization  to  pXOl  genes,  including 
all  three  anthrax  toxin  genes  (Fig.  6,  which  is  published  as 
supporting  information  on  the  PNAS  web  site). 

Taken  collectively,  these  data  provided  multiple  lines  of 
evidence  to  suggest  that  the  G9241  isolate,  responsible  for  an 
anthrax-like  clinical  presentation,  is  indeed  a  B.  cereus  isolate 
containing  pXOl  genes,  including  the  genes  encoding  the  an¬ 
thrax  toxins. 

Genome  Sequence  Analysis.  To  obtain  a  more  in-depth  under¬ 
standing  of  this  unusual  B.  cereus  clinical  isolate,  whole-genome 


B.  cereus  B.  thuringiensis 

ATCC 14579  tolworthii  B-  cereus  B.  cereus 


Fig.  2.  Unrooted,  neighbor-joining  tree  derived  from  multiple-locus  se¬ 
quencetyping  of  6.  cereus  G9241  and  other  Bacillus  spp.  (http://pubmlst.org/ 
bcereus). 
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shotgun  sequencing  was  performed.  One  feature  of  the  B.  cereus 
G9241  genome  is  its  mosaic  structure.  This  structure  could  be 
explained  by  the  presence  of  a  great  number  of  known  and  novel 
putative  mobile  genetic  elements  (insertion  sequence  elements, 
phages,  transposases,  integrases,  and  recombinases).  Fig.  5A 
shows  the  degree  of  homology  of  each  of  the  B.  cereus  G9241 
predicted  proteins  to  either  B.  anthracis  Ames,  B.  cereus 
ATCC10987,  or  B.  cereus  ATCC14579.  The  proteome  of  B. 
cereus  G9241  shows  a  greater  degree  of  similarity  to  the  core 
proteome  of  the  complete  genome  of  B.  cereus  ATCC10987  than 
to  B.  cereus  ATCC14579  and,  to  a  lesser  extent,  B.  anthracis 
Ames,  also  shown  by  blast  score  ratio  analysis  (Fig.  5 B). 
Overall,  60%  of  the  G9241  proteome  is  more  similar  to  B.  cereus 
strains  than  to  B.  anthracis,  in  full  agreement  with  the  initial 
classification  of  G9241  as  a  B.  cereus  strain. 

Comparative  analysis  of  the  chromosome  of  B.  cereus  G9241 
with  those  of  B.  cereus  ATCC10987,  B.  cereus  ATCC14579,  and 
B.  anthracis  Ames  revealed  a  set  of  specific  characteristics  that 
can  be  linked  to  phenotypic  observations.  PlcR  is  a  positive 
regulator  of  many  secreted  and  virulence-associated  genes  in  B. 
cereus  strains,  including  protease,  phospholipases,  hemolysins, 
and  motility  genes  (24).  In  common  with  other  B.  cereus  stains, 
G9241  does  not  have  the  nonsense  frameshift  found  in  the  plcR 
gene  of  B.  anthracis',  hence,  it  could  produce  a  fully  functional 
PlcR  regulator,  which  may  explain  many  of  its  phenotypic 
characteristics  (Fig.  7,  which  is  published  as  supporting  infor¬ 
mation  on  the  PNAS  web  site). 

B.  cereus  G9241  is  motile  (Table  1),  a  phenotype  that  could  be 
explained  by  the  presence  of  a  complete  flagella  biosynthetic 
cluster.  Although  B.  anthracis  is  not  motile,  it  is  common  for  B. 
cereus  and  other  motile  Bacillus  species  to  be  flagellated.  Even 
though  B.  anthracis  contains  most  of  this  system,  four  essential 
proteins  are  truncated  by  frameshifts  (10),  all  of  which  are 
full-length  in  G9241.  Two  PlcR  palindromic  binding  sites  are 
located  124  bp  in  front  of  a  flagella  component  subunit  ( flgB ) 
and  a  conserved  hypothetical  protein  that  is  part  of  the  same 
cluster.  Interestingly,  this  isolate  contains  five  flagellin  subunits 
(only  two  are  found  in  B.  cereus  ATCC10987,  four  in  B.  cereus 
ATCC14579,  and  one  in  B.  anthracis  Ames),  four  of  which  are 
clustered  and  have  low  similarity  to  either  B.  anthracis  or  B. 
cereus  flagellins.  However,  the  fifth  is  most  similar  to  a  flagellin 
subunit  that  appears  to  be  only  found  in  B.  cereus  (Fig.  8,  which 
is  published  as  supporting  information  on  the  PNAS  web  site). 

The  genome  of  G9241  encodes  for  four  putative  hemolysins 
(hemolysin  A,  BL,  II,  and  III).  Hemolysin  BL  appears  to  be 
PlcR-regulated  as  in  other  B.  cereus  strains,  which  could  explain 
the  observed  hemolytic  activity  of  the  G9241  isolate  (Table  1). 
In  addition,  G9241  phospholipase  C,  which  is  responsible  for 
lecithin  degradation,  has  all  its  components  under  the  control  of 
PlcR  as  do  other  B.  cereus  strains. 

An  unexpected  finding  in  the  G9241  genome  was  the  presence 
of  a  191,110-bp  circular  plasmid  with  high  similarity  and  synteny 
to  B.  anthracis  pXOl  (Fig.  1).  This  plasmid  was  responsible  for 
the  high  level  of  hybridization  to  pXOl  genes  observed  by  CGH. 
We  will  refer  to  this  plasmid  as  pBCXOl  to  differentiate  it  from 
B.  anthracis  pXOl.  pBCXOl  encodes  for  the  entire  anthrax  toxin 
biosynthetic  complex,  PA  (99.7%  amino  acid  identity),  lethal 
toxin  (99%  amino  acid  identity),  and  edema  factor  (96%  amino 
acid  identity)  and  its  associated  regulatory  proteins  AtxA  (100% 
amino  acid  identity)  and  PagR  (98.6%  amino  acid  identity), 
confirming  the  CGH  and  PCR  results.  The  presence  of  AtxA  and 
a  full  PlcR  regulon  raises  some  interesting  questions,  because  it 
was  suggested  previously  that  in  B.  anthracis  these  two  regulators 
were  incompatible  (25)  and  may  have  led  to  the  selection  for  the 
plcR  nonsense  mutation  in  B.  anthracis.  The  B.  cereus  G9241 
pBCXOl-encoded  pagA  gene  represents  a  distinctive  PA  se¬ 
quence  type  by  using  the  typing  scheme  of  Price  et  al.  (26).  Two 
point  mutations  are  found  at  positions  1,999  and  2,672,  resulting 


Fig.  3.  B.  cereus  G9241  capsule.  India  ink  stain  of  6.  cereus  G9241  grown  in 
ambient  atmosphere  (A)  and  B.  anthracis  grown  in  5%  C02  (6).  (Scale  bar, 
««10  |ixm.) 

in  serine  to  proline  and  isoleucine  to  serine  amino  acid  changes, 
respectively  (Table  4,  which  is  published  as  supporting  informa¬ 
tion  on  the  PNAS  web  site).  With  these  two  mutations,  the 
identified  sequence  is  most  similar  to  PA  genotype  V,  which  has 
been  primarily  associated  with  the  western  North  America 
diversity  group  of  B.  anthracis  (26,  27).  Expression  of  the  PA  by 
B.  cereus  G9241  was  confirmed  in  the  supernatant  of  a  stationary 
phase  culture  (109  cells  per  ml)  by  ELISA  (see  Materials  and 
Methods).  The  amount  of  PA  expressed  (4.3  jug/ml)  was  found 
to  be  equivalent  to  that  seen  in  B.  anthracis  Ames  (3.4  pg/ml) 
and  B.  anthracis  Sterne  (2.7  p,g/ml).  Importantly,  the  negative 
controls  B.  anthracis  Pasteur  and  B.  cereus  ATCC10987  did  not 
express  any  detectable  amount  of  PA. 

By  using  a  two-component  DFA  assay  (16),  B.  cereus  G9241 
was  tested  (Table  1)  for  cell  envelope  similarities  to  B.  anthracis. 
G9241  produced  a  positive  reaction  with  mAb  specific  to  the  B. 
anthracis  cell-wall  antigen  but  not  with  Ab  against  the  B. 
anthracis  capsule.  Because  it  is  not  uncommon  for  some  clinical 
isolates  of  B.  cereus  to  react  with  the  cell-wall-specific  Ab,  the 
positive  reaction  with  this  Ab  was  not  unusual  (16).  Although  the 
isolate  did  not  react  with  the  B.  anthracis-specific  capsule  Ab,  it 
did  produce  a  capsule  that  was  expressed  during  growth  on  sheep 
blood  agar  plates  or  in  horse  serum  (Fig.  3)  but  did  not  require 
CO2  induction,  as  is  the  case  for  B.  anthracis.  However,  a  thicker 
capsule  was  observed  during  growth  with  elevated  CO2  (M. 
Chute,  personal  communication).  A  novel  polysaccharide  cap¬ 
sule  biosynthetic  gene  cluster  has  been  found  on  a  unique  and 
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Fig.  4.  Survival  of  A/J  mice  i.p.  challenged  with  B.  cereus  G9241  (green),  6. 
anthracis Sterne  (yellow),  and  B.  cereus  ATCC1 0987  (red).  High-spore  dose  (1  X 
106)  and  low-spore  dose  (1  X  1 04)  are  represented  by  triangles  and  circles, 
respectively.  The  experiment  was  monitored  for  14  days,  after  which  the  mice 
inoculated  with  6.  cereus  ATCC1 0987  were  still  alive. 


previously  unidentified  218,094-bp  circular  plasmid  in  the  ge¬ 
nome  of  G9241.  pBC218  encodes  188  putative  ORFs,  among 
which  we  identified  gene  products  with  similarities  to  B.  anthracis 
pXOl-encoded  toxin  regulator  AtxA  (78%  amino  acid  identity), 
PA  (60%  amino  acid  identity),  and  the  lethal  factor  (36%  amino 
acid  identity)  but  not  the  edema  factor  or  the  B.  anthracis 
pX02-encoded  poly-y-D-glutamic  acid  capsule  biosynthetic 
genes.  The  pBC218-encoded  PA  peptide  is  identical  to  B. 
anthracis  at  all  10  dominant  negative  amino  acid  residues  iden¬ 
tified  by  Mourez  et  al.  (28).  Another  33  residues  were  identified 
that  resulted  in  decreased  activity  when  substituted  by  cysteine, 
of  which  27  are  identical  in  pBC218-encoded  PA.  The  remaining 
six  residues  are  conservative  substitutions  (28).  This  evidence 
suggests  that  the  pBC218-encoded  PA  could  be  functional; 
however,  this  must  be  verified  experimentally.  Because  no 
pX02-like  plasmid  was  identified  in  the  genome,  the  encapsu¬ 
lation  of  G9241  observed  in  the  laboratory  (Table  1)  may  be 
encoded  by  the  pBC218  capsule  biosynthetic  cluster.  Analysis  of 
this  cluster  indicates  that  it  encodes  for  a  polysaccharide  capsule 
(Fig.  9,  which  is  published  as  supporting  information  on  the 
PNAS  web  site).  In  addition,  our  detailed  genomic  examination 
has  not  revealed  any  other  capsule  biosynthetic  clusters  on  either 
the  chromosome  or  plasmids. 

Lethality  of  B.  cereus  G9241  in  Mice.  These  genotypic  characteristics 
of  the  G9241  B.  cereus  strain  supported  by  the  severe  clinical 
presentation  are  consistent  with  its  properties  as  a  B.  anthracis- 
like  pathogen.  To  assess  the  mammalian  virulence  of  B.  cereus 
G9241,  A/J  mice  were  challenged  by  i.p.  injection  (21)  with 
mature  spores  of  B.  anthracis  Sterne,  B.  cereus  G9241,  and  B. 
cereus  ATCC10987  (12).  In  these  preliminary  studies,  G9241 
produced  100%  lethality,  consistent  with  the  positive  control  B. 
anthracis  Sterne,  although  lethality  was  achieved  faster  when  the 
high-spore  dose  was  inoculated  (f  test,  P  =  0.081)  (Fig.  4).  B. 
cereus  G9241  colonies  were  recovered  from  the  heart  tissue,  and 
the  blood  was  removed  from  the  animals  postmortem  and  plated 
on  sheep  blood  agar  incubated  at  37°C  for  24  h.  Bacterial 
colonies  were  not  recovered  from  any  of  the  negative  control 
mice  injected  with  either  PBS  or  with  the  nonpathogenic  B. 
cereus  ATCC10987. 


Additional  Genome  Features.  Several  other  interesting  features  of 
this  strain  included  a  cryptic  phage  of  29,886  bp  (pBClin29,  Table 
2)  that  encodes  phage-like  proteins  and  a  plasmid  replicon 
similar  to  B.  anthracis  pX02.  In  addition  and  particularly 
relevant  to  public  health  and  patient  treatment,  G9241  contains 
genes  that,  if  functional,  should  provide  the  isolate  with  resis¬ 
tance  to  j3-lactam,  chloramphenicol,  and  macrolide  antimicro¬ 
bial  agents  (Table  4).  A  PlcR  palindromic  binding  site  is  located 
48  bp  in  front  of  a  metallo-j3-lactamase  (ORF04823),  which 
could  explain  the  penicillin-resistant  phenotype  observed  in  the 
laboratory.  This  binding  site  is  also  present  in  the  other  B.  cereus 
genome  and  in  B.  anthracis  Ames.  However,  the  truncated  PlcR 
is  proposed  to  explain  the  penicillin-susceptible  phenotype 
observed  in  the  latter. 

Discussion 

This  article  describes  the  detection  of  naturally  occurring  anthrax 
toxin  genes  in  a  species  other  than  B.  anthracis  and  associated 
with  an  inhalation  anthrax-like  illness.  Although  identification  of 
B.  cereus  as  the  etiologic  agent  of  pneumonia  is  not  common,  it 
has  been  reported  and  is  generally,  but  not  exclusively,  associated 
with  significant  risk  factors.  Miller  et  al.  (22)  reported  on  two 
fatal  cases  of  fulminating  bacteremia  and  pneumonia  that  oc¬ 
curred  in  1995  and  1996  with  similarity  in  their  clinical  presen¬ 
tation  to  the  present  study,  but  those  B.  cereus  strains  did  not 
contain  pagAAike  genes.  Recently,  two  cases  of  fatal  respiratory 
illness  in  the  United  States  were  also  attributed  to  B.  cereus 
isolates,  and  preliminary  tests  indicated  the  presence  of  a  B. 
anthracis -like.  pagA  gene  in  these  isolates  (A.R.H.  and  T.P., 
unpublished  data).  Further  characterization  of  these  and  other 
B.  cereus  isolates  from  cases  of  severe  pneumonia,  for  which  the 
Centers  for  Disease  Control  and  Prevention  is  in  the  process  of 
setting  up  a  prospective  surveillance,  could  potentially  reveal  the 
true  burden  of  disease  caused  by  isolates  harboring  what  are 
traditionally  considered  B.  anthracis  virulence  plasmids  and/or 
their  genes.  The  presence  of  the  pXOl  and  pX02  plasmid- 
associated  genes,  however,  is  not  normally  investigated  in  an 
isolate  unless  it  is  suspected  to  be  B.  anthracis.  In  addition,  the 
isolation  of  B.  cereus  in  a  clinical  laboratory  is  frequently 
considered  contamination,  resulting  in  it  being  discarded  rather 
than  to  follow  up  with  further  characterization  and  reporting  (5). 
Consequently,  clinical  isolates  of  species  other  than  B.  anthracis 
containing  pXOl  or  its  genes  may  have  been  missed. 

The  rapidity  with  which  whole-genome  sequencing  data  can 
currently  be  obtained  provides  a  new  approach  to  the  public 
health  challenges  associated  with  emerging  diseases  (29).  The 
high-coverage  draft  genome  data  on  G9241  not  only  clarified  the 
relationship  of  this  isolate  to  other  members  of  the  B.  cereus 
group  but  also  provided  direct  insight  into  the  correlation  of 
phenotypic  characteristics  and  the  genetic  basis  for  their  pre¬ 
sentation.  Recent  studies  have  shown  that  genes  that  have 
traditionally  been  considered  restricted  to  B.  cereus  or  B.  an¬ 
thracis  are  actually  more  widespread  (10,  30-32).  Although 
genomic  sequences  with  homology  to  B.  anthracis  pXOl  and 
pX02  sequences  have  been  detected  in  other  Bacillus  spp.  (12), 
this  work  reports  on  the  toxin  genes  and  an  entire  pXOl-like 
plasmid  found  outside  of  a  B.  anthracis  isolate. 

Identification  of  Bacillus  spp.  containing  the  B.  anthracis 
virulence  genes  and  able  to  cause  severe  respiratory  illness  that 
could  potentially  be  confused  with  inhalation  anthrax  adds 
substantially  to  the  complexity  of  clinical  and  laboratory  diag¬ 
nosis  in  general  and  particularly  during  a  potential  bioterrorism 
event.  Given  that  pXOl-  and  pX02-based  genes  and/or  their 
encoded  products  are  the  most  frequently  used  targets  to  detect 
and  identify  B.  anthracis,  false  positive  results  could  occur  and 
result  in  unnecessary  alarm.  On  the  other  hand,  it  may  not  be 
appropriate  to  consider  B.  anthracis,  as  currently  defined,  as  the 
only  species  capable  of  causing  inhalation  anthrax-like  disease. 
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Although  this  isolate  caused  a  respiratory  illness,  it  will  also  be 
interesting  to  study  the  potential  of  this  isolate  to  cause  gastro¬ 
intestinal  or  cutaneous  infections  similar  to  B.  anthracis.  Limiting 
screening  for  B.  anthracis- specific  phenotypes  rules  out  isolates 
that  are  hemolytic  and/or  motile,  such  as  the  B.  cereus  described 
in  this  study.  A  more  sophisticated  system  to  prospectively 
monitor  for  the  presence  of  B.  anthracis- specific  virulence 
plasmids  and  plasmid-associated  genes  in  Bacillus  spp.  other 
than  B.  anthracis  may  have  to  be  developed.  Vigilance  and  close 
communication  between  clinicians  treating  patients  and  labora¬ 
tory  scientists  is  needed  to  assure  that  the  true  causative  agents 
are  rapidly  isolated  and  identified. 

The  discovery  of  a  pathogen  such  as  B.  cereus  G9241  also 
emphasizes  the  challenges  of  creating  lists  of  potential  biothreat 
agents.  Although  it  is  not  practical  to  include  every  pathogenic 
species  on  such  lists,  it  should  be  noted  that  depending  on  the 
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